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This document presents a ray-optics analysis of the response of a two-element radio interferometer 
to the presence of a plane gravitational wave. A general expression for the differential phase observed 
between the two receiving stations as a result of an arbitrary gravitational wave is determined, as 
well as the specific responses to monochromatic and black hole ring-down waveforms. Finally, 
the possibility of gravitational wave detection via this mechanism is discussed in the context of 
interferometer noise. 
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I. INTRODUCTION 

A gravitational wave is a perturbation to the geometry 
of space that propagates through space at the speed of 
light. It's effect is to stretch and shrink the distances be- 
tween objects in the directions transverse to its direction 
of travel. The distortion is volume-preserving, so while it 
is shrinking distances in one direction, it stretches them 
in another. 

The stretching and shrinking effect is not a co-ordinate 
artifact. A measuring stick used to monitor the distance 
between two objects freely-falling in the presence of a 
gravitational wave will measure one separation distance 
at one time and a different separation later. The stretch- 
ing and shrinking is also linear: given the change in dis- 
tance between two objects separated by a distance L, the 
change in distance between two objects separated by 2L 
will be doubled. The amplitude of a gravitational wave is 
expressed as a dimensionless quantity giving the change 
in distance per unit distance — a shear strain in space. 

Any apparatus capable of accurately measuring the 
distance between two objects, or the change in the dis- 
tance between two objects, as a function of time is, in 
principle, capable of detecting a gravitational wave. The 
question of the effectiveness of the apparatus in accom- 
plishing this task is one of sensitivity and signal-to-noise 
ratio. One design being actively pursued by several re- 
search groups today is based on a Michelson-type laser 
interferometer. In this design, a laser beam is split and 
sent down the two orthogonal arms of the apparatus, 
and then at the end of each is reflected off a suspended 
mass back toward the splitter where the beams are re- 
combined 

0, S II II 13 If the length of one or the 
other arms changes, the result will be a change in the 
differential phase between the returning laser beams and 
so a passing gravitational wave will appear as a shift 
in the interference pattern at the output of the inter- 
ferometer. A variation on this technique that is also 
being pursued by several research groups is a detector 
consisting of an Earth-based transmitter/receiver and a 
spacecraft carrying a radio or optical transponder. There 




FIG. 1: A gravitational wave moving across the lines of sight 
in the plane of a two-station interferometer. One of the gravi- 
tational wave's polarization axes is aligned (ignoring parallax) 
with the lines of sight. 



are currently proposals to launch specially-designed ve- 
hicles dedicated to this task usingoptical links from the 
ground to the vehicle and back [6(. It is also possible 
to use existing missions. Deep-space vehicles generally 
carry radio transponders that are used as aids to naviga- 
tion. By maintaining a coherent up- link/down- link chan- 
nel at high frequency stability with appropriate modula- 
tions one can measure the range and range-rate to the 
spacecraft with high precision and thereby search for the 
assage of a gravitational wave through the solar system 
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In this document, another type of gravitational wave 
detection apparatus similar to this last technique is intro- 
duced and investigated. Figure ^ shows a radio interfer- 
ometer arrangement with a distant radio source and two 
receiving elements. It also shows a gravitational wave 
travelling across the lines of sight to the radio source in 
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the plane defined by the source and the two interferome- 
ter elements (to be referred to, hereinafter, as the plane of 
the interferometer). For the purpose of illustration, one 
of the gravitational wave's polarization axes is aligned 
(ignoring parallax) with the lines of sight. Because the 
gravitational wave travels at finite speed, any perturba- 
tion to distances transverse to its direction of travel will 
influence the signal arriving at receiver 1 in advance of its 
influence on the signal arriving at receiver 2. The pres- 
ence of a gravitational wave will, therefore, manifest itself 
as a time-varying difference between the arrival times of 
radio signals at the receivers — a time-varying differen- 
tial phase between the stations. Such a detector is similar 
in operation to the single arm, spacecraft-based, coherent 
up-link/down-link type of detector. In this case, however, 
a second arm is introduced to act as a phase reference for 
the first which permits the use of a broad-band radio sig- 
nal. What follows is an analysis of the magnitude of the 
differential phase between the stations when a gravita- 
tional wave is present, and a preliminary assessment of 
the signal-to-noise ratio in this type of detector. 

We will consider a radio signal emanating from the ra- 
dio source toward each of the receiving stations. Each 
radio signal is treated as a null geodesic in space-time, 
originating on the radio source's world line. The inter- 
sections of the geodesies with the world lines of the re- 
ceivers define events in space-time giving the co-ordinate 
times of arrival of the signal at each receiving element. 
By determining the difference in the times of emission of 
the signals being received simultaneously at the two re- 
ceiving elements, and multiplying this time difference by 
the angular frequency of the radio emission one obtains 
the differential phase between the two receiving elements. 
The calculation will first be done for arbitrary receiver 
locations and arbitrary gravitational waveforms, then a 
simplified interferometer arrangement will be introduced 
in order to more easily discuss the response to specific 
gravitational waves. Throughout the calculations, the 
units for length will be chosen so as to make the numer- 
ical value of the speed of light, c, equal to 1. 



II. NULL GEODESICS IN A PLANE 
GRAVITATIONAL WAVE SPACE-TIME 

For simplicity, we assume a Minkowski background 
for the geometry of space-time. Rather than solving 
the geodesic equation for fixed radio source and receiver 
positions in an arbitrarily-oriented gravitational wave, 
we shall choose to fix the orientation of the gravita- 
tional wave and adjust the geodesic's end-points. The 
co-ordinate system is defined as follows: the radio source 
locates the origin of the spatial co-ordinates, the gravita- 
tional wave vector sets the negative- z direction and the 
two polarization axes of the gravitational wave set the x 
and y directions. The origin of the time co-ordinate is ar- 
bitrary. Space-time geometry is described by the metric 



tensor 

SV = V + V" W 
where the background Minkowski metric tensor is 

V, = diag(-1, 1,1,1), (2) 

and the gravitational wave is described by the tensor h^ v 
given by 



(3) 



e+, is the gravitational wave's polarization tensor which, 
in this case, is 



diag(0, 1,-1,0), 



(4) 



and h(z + t) is the gravitational wave's amplitude and 
describes a waveform propagating at the speed of light 
in the negative z direction. h(z + t) carries units of strain 
- change in length per unit length making it dimen- 
sionless. Since the radio source defines the origin of the 
spatial co-ordinates, the calculations are being done in its 
rest frame so co-ordinate time corresponds to the proper 
time measured by a clock carried by the radio source. 
This may or may not agree with the proper times mea- 
sured at the receiving elements. 

We assume the amplitude of the gravitational wave is 
small, | h | <C 1, and carry out calculations to first order in 
h. Altogether, the line element describing the geometry 
of space-time is 

ds 2 = - dt 2 + [1 + h(z + t)} dx 2 + 

[l-h(z + t)]dy 2 + dz 2 . (5) 

The only non-zero Christoffel symbols for this line ele- 
ment are 



t yy = ^h'(z + t) 



r 

pa; 



= - TZ yy = -l h '( z + t ) 
tx =T\ t = l -h'{z + t) 

zx = T* xz = \h\z + t) 



ty 



vy yt = -\h\z + t) 



ry zy = ry vz = --h'(z + t), 



where h' is the derivative of h with respect to its ar- 
gument: h'(z + t) = d ( z d ft ) fe(^ + t). A geodesic in the 
space-time is a solution of the geodesic equation, 



d 2 x» 
dA 2 



dx a dx 

Q/3 ~dT~dT 



o. 



(6) 



and is described by four functions of an affine parame- 
ter: x(A), y(\), z(X), and i(A), where A parameterizes 
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events along the the geodesic. Using the Christoffel sym- 
bols from above in the geodesic equation leads to the 
following four coupled, second-order, non-linear, differ- 
ential equations for x(X), y(X), z(X), and t(X) along a 
geodesic in this space-time: 



dH 
dX 2 



1 



h'(z + t) 



dx 
dA 



dy 
dX 



= 0, 



dA 2 "^ 



d2 y h '( 



dz 
dA 
dz 
dA 



d* 
dA 
dt 
dA 



d 2 z 1 



da; 
dA 



dy 
dA 



(7a) 

(7b) 
(7c) 

(7d) 



A unique solution to these equations requires eight 
boundary /initial conditions, seven of which are provided 
by our knowledge of the end points of the geodesies: we 
know the spatial co-ordinates of the radio source and the 
spatial co-ordinates of the receiving element and also ei- 
ther the time of emission or the time of reception. If 
we choose A = to correspond to the emission event 
and A = 1 the reception event on the radio signal's null 
geodesic, then the boundary conditions are 



t(X 


= 0) 


— to, 


t(X = 


l)=t u 


(8a) 


x(X 


= 0) 


= o, 


x(X = 


1) = x i: 


(8b) 


y(A 


= 0) 


= 0, 


2/(A = 


1) = Vi, 


(8c) 


z(A 


= 0) 


= 0, 


z(X = 


1) = Zi, 


(8d) 



where Xi, yi, Zi, and U are the space-time co-ordinates 
of the reception event at the i th receiving element, and 
to is the co-ordinate time of emission of the radio signal. 
The subscripts on the reception event's co-ordinates will 
be used later to enumerate them when multiple receiving 
elements are considered. We have introduced both to 
and U but only one of the two is known. The other is 
to be found from the solution to the geodesic equation. 
The eighth constraint on the solution will be provided by 
imposing the requirement that the geodesic be null. 
Adding {Taj to ||7dJ| gives 



dA 2 



(z + t) = 



(9) 



which requires that whatever form z(X) and t(A) take, 
their su m must be a linear function of A. We introduce 
r, = \J x 2 + y 2 + z 2 to be the distance from the radio 
source to the i th receiving element in the absence of the 
gravitational wave. This provides the time-of- flight, ti — 
to, for the radio signal from source to receiver in the 
absence of the gravitational wave. The presence of a 
gravitational wave alters this time-of-flight to the new 
value a\Ti where a\ is an unknown dimcnsionless quantity 
to be found from the solution to the geodesic equation. 
The co-ordinate time of emission, to , of the radio signal 



can be expressed in terms of the the co-ordinate time 
of reception and the time-of-flight in the presence of a 
gravitational wave as 



t = t l - a\T{. 



(10) 



Using this and the boundary conditions on z(A) and i(A), 
we can write the solution to I© as 

, ( A) + t(A)^{^+ a ^+;°' (id 

where it is expressed in terms of either t or ti (depend- 
ing on which of the two is known), and the unknown a\ 
which appears here as an integration constant. Having 
the solution in l(TT|) succeeds in decoupling each of the 
equations for x(A) and y{X) in Q from the other three 
since it is only z + t and its derivatives that appear in 
those two equations. 
Introducing 

r(A)=z(A)+*(A), (12) 
we can write h'(z + t) as 

^ + *) = ^Mr)=(^)" 1 ^[r(A)]. 

From Qllfl. dr / dA is a constant (independent of A) equal 
to 



dr dz dt 

-x = Jx + -x =Zl + ain - 

The x equation in l|7b|l then becomes 
d 2 x dx d 



(13) 



dA 2 dA dA 



h(r) = 



which can be integrated once to get 

= Xi [l + a 2 ]e-' l M a Xi [l + a 2 ][l - h(r)}, 
and then integrated again to get 



x(X) = Xj 



1 + a 2 



A- 



^[T(A')]dA' + a3 , (14a) 



where a 2 and 03 are constants of integration. The bound- 
ary conditions on a: (A) require 



a 2 = f h[T(X')]dX', 
Jo 



a 3 = 0. 

Similarly, the y equation in l)7c|l becomes 
d 2 y dy d 



(14b) 
(14c) 



dA 2 



dA dA 



Kt) = 0, 
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and by comparison with the x equation and its boundary 
conditions, it can be verified that the solution for y(A), 
satisfying its boundary conditions, is 



2/(A) = Vi 



1 - a 2 



X + / h[r(X')]dX' -a 3 
Jo 



(15) 



where ai and 0,3 are the same integration constants as in 
the solution for x(X). 

For z(X) in |7dJ, we need the derivatives of x(X) and 
y(X) but we only need them to th order in h because 
higher order terms will result in quadratic and higher or- 
der terms when multiplied by h! in the differential equa- 
tion. These derivatives are 



Vi 



dx dy 
~dX~ Xl ' dA 

With these, the z equation in l|7d(l becomes 



d 2 z 



y, 



dX 2 2 V dr/ dA j dX 
which can be integrated twice to give 



h{r) = 



z(X) = z l 



A+- 



1 x 2 - y 2 



2(dT/dX) Zi J 



/i[r(A')] dA' - a 4 A - a 5 



(16a) 

where 04 and 05 are integration constants. The boundary 
conditions on z (A) require 



1 x 



2(dr/dA)z i y 
a 5 = 0. 



h[T(X')]dX', 



(16b) 
(16c) 



The eighth and final constraint on the solution is ob- 
tained from the requirement that a space-time interval 
along the radio signal's geodesic be null. Setting ds 2 = 
in jnj yields 



dt 
dA 



[l + h{z + t)\ 



dx 
dA 



[l-h{z + t)\ 



dy 
dA 



dz 
dA 



and using (|T2l for dt/ dA, this becomes 
2 



dr 
dA 



dr dz 
dA dA 



d.?' 



| 1 + „ MI| _ + [i-*W]UJ 



dy 



(17) 



Now, to first order in h 

dx 

dA = 
dy _ 
dA 



1 - H(t) 
1 + H(t) 



dz 
dA 



1 + \4vm> HiT) 



where 



H(t) = h(r) 



h[T(X')]dX'. 



With these, the constraint in (|17H becomes 



dr \ „ dr , , 9 _ TT 
dA - 2 dA 2i "^-^ 



x 2 [l + h(r) - 2H] + y 2 [l - h(r) + 2H] 



or 



2 2 ci 



h[T(X')]dX'. 



Using (|13|l for dr/dX and recalling that there are two 
expressions for t depending on which end-point of the 
geodesic is known, this becomes 

a\ = \+ 1 2 Vl / h[T(X')]dX' 
r i Jo 

1 + Jo h[{ Zi + a in )X' + t Q ] dA', 

2 2. 2 1 

1 + ^r 1 - J h[(zi + airi)X' + U - aifj] dA'. 

This constitutes a transcendental equation for the inte- 
gration constant a\ — the last unknown in the solution 
for the geodesic. This equation may be solved by recog- 
nizing that the right-hand side is equal to 1 plus a quan- 
tity proportional to h, and so requiring the left-hand side 
to also be equal to 1 plus a quantity proportional to h. 
Consequently it is convenient to write 



1 

ai = 1 + -a 6 



(18a) 



where must be proportional to h. Making this substi- 
tution and discarding terms higher than first order in h 
yields 1 



t^ltfhKzi + rjX' + to] dA', 
El^L^h[(z l + r % )X' + t l ~r l ]dX'. 



A change of variables in the integrals from A' to r turns 
these expressions for ae into 



Xi ~ Vi Mr) dr, 



„ J (Zi+ri)rf Jt 

a 6 — \ „2 .2 



(18b) 



(zi-\-ri)r? Jti — ri 



1 The a§s inside the integrand on the right-hand side are dis- 
carded. To do this carefully, one expands the integrand in a 
Taylor series in powers of ag. Since ag must be proportional to 
h, only the zeroth-order term in that series need be retained as 
all higher order terms are quadratic or higher in h. 
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Finally, using the expressions for ae in a\ and substi- 
tuting the results into (|10H we obtain 



U — in 



2(zi + n)n 



h(r)dT, (19) 



and 



to = U — n — 



Vi 



2(zi + r i )r l 



U+Zi 



h(r) dr. 



(20) 



These give, respectively, the co-ordinate time of reception 
as a function of the receiver position and the co-ordinate 
time of emission, (|19|l , and the co-ordinate time of emis- 
sion as a function of the receiver position and the co- 
ordinate time of reception, l|20() . From these we see that 
the time-of-flight of the radio signal is given by plus 
a term expressing the perturbation of the radio wave's 
travel time due to its passage through the gravitational 
wave en route to the point of reception. 








a 
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III. THE DIFFERENTIAL PHASE IN A 
TWO-ELEMENT RADIO INTERFEROMETER 

Consider an interferometer system consisting of two re- 
ceiving elements and a data acquisition system at each 
element which is capable of establishing co-ordinate time 
and recording the radio signal being received as a func- 
tion of this parameter. In general, the signals being 
received at the same co-ordinate time at each element, 
ti = t 2 = t, will have been emitted from the radio source 
at different co-ordinate times. The difference in the times 
of emission can be obtained from (|20() and is 



Ai = r 2 - n- 

1 

2 



2 2 
X 2 ~ V2 



(z 2 + r 2 )r 2 J t - 



rt+z 2 

I h{r) dr- 

Jt-r 2 



2 2 

x{ - yf 



(zi + n)ri J t _ ri 



t+zi 



h(r) dT 



r 2 — T\ is the flat-space geometric delay for the interfer- 
ometer and the remainder of the expression represents 
an additional contribution to the differential emission 
time arising from the presence of a gravitational wave. If 
the radio receiving elements are tuned to an angular fre- 
quency of u> T (in the rest frame of the radio source), then 
this differential emission time translates into a differen- 
tial phase of the received signals given by A</> = oj T Ato, 
or 



A0(t) = uj,(t 2 - n)+ 



~2 



vi 



{z 2 + r 2 )r 2 



t + Z 2 



h(r) dr- 



f— r 2 



(zi + r x )ri 



t+Zl 



h{r) dr 



(21) 



FIG. 2: The orientation of the interferometer with respect to 
the gravitational wave's wave-vector and polarization. 



The differential phase in l|21|) is quite general. In prin- 
ciple, for example, the expression continues to hold for 
radio receiving elements that are in motion with respect 
to the radio source and will correctly account for effects 
such as the Doppler shifting of the radio signals. The 
expression, however, is arrived at in the rest frame of 
the radio source. This means, for example, that the co- 
ordinate time, ti, at each receiving element is the proper 
time measured by a clock carried by the radio source, not 
by a clock carried by that receiving element. In practise 
it is local proper time that is available at each receiving 
element. The implementation of this gravitational wave 
observation technique will require knowledge of the trans- 
formation from local proper time to co-ordinate time. For 
a discussion of this transformation, see 0. 

We now consider a special case for the interferome- 
ter configuration consisting of two receivers at rest with 
respect to and equally distant from the radio source, 
r\ = r 2 = r, and separated from each other by a distance 
21 where l/r -C 1. The geometry of the arrangement 
is depicted in Figure [2 We introduce three angles to 
describe the interferometer's orientation with respect to 
the gravitational wave: two to describe the direction from 
which the gravitational wave is approaching and a third 
to describe the orientation of its polarization axes. The 
direction from which the gravitational wave is approach- 
ing is defined by (9,(f>), the polar angle and azimuthal 
angle of a spherical polar co-ordinate system. The polar 
axis, 9 = 0, is directed toward the radio source, bisecting 
the interferometer's baseline. Also bisecting the inter- 
ferometer's baseline and perpendicular to the polar axis 
is the origin of the <fi co-ordinate. The orientation of 
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the gravitational wave's polarization tensor is given by a 
which is measured in the plane orthogonal to the gravita- 
tional wave vector at the point (0, (f>), positively from the 
meridian of constant <$> to the polarization axis as shown 
in Figure |2 The choice of which polarization axis is used 



to define the angle a is arbitrary since the two possible 
choices simply differ from one another by a rotation of 
7r rad in the detector's output. 

In the co-ordinate system used to obtain the differen- 
tial phase in it2*T]) . the receiver positions are given by 



cos a sin a ' 
— sin a cos a 
1, 



' cos — sin 
1 
sin cos 




sin cos a ± - 
— sin 9 sin a =F 

— cos 



(sin 4> cos 9 cos a 
- (sin <p cos " 
± l - 



sm a 
sin d> sin i 



cos (j) sin a) 
- cos 4> cos a) 



The rotation matrices are those that will rotate the grav- 
itational wave vector in Figure [3 back to the z axis while 
orienting its polarization directions along the x and y 
axes. From this, 



2 2 

x -y 



r z sin z 9 cos 2a\ 1 ± 2-- 



l cos 9 sin < 



and 



[z + r)r 



z + r = r(l — cos0)(^ 



= (1 + cos 9) cos 2c 



1± 



r sin 9 cos 2a 
I sin 6 sin 



r 1 — cos ( 



I / cosflsincA sinc6sin< 

1 ± -12 

r V sin 9 cos 2a 1 — cos 6 



Using this expression in (|21H and considering terms in 
h- to be second order and negligible, the interferometer's 
differential phase reduces to 



Ac6(t) = y(l + cos6»)cos2a 



t—r cos 9-\-l sin sin (ft 



h{r) dr. 



t—r cos # — £ sin # sin c 



Redefining the origin of the time co-ordinate so that 
t — rcos9 — > t, the differential phase for the simplified 
interferometer configuration is finally 



A<j>(t) 



2 



/"' 4 

(1 + cos 9) cos 2a I 

Jt-i 



t+l sin 8 sin < 



/i(r)dr. (22) 



sin V sin < 



We can also express the differential phase in terms 
of the gravitational wave's frequency-domain represen- 
tation. Introducing the Fourier transform of the gravita- 
tional wave, 



h(r) 



h{oj) = 



+ 00 



h(tu)e iLUT dw, 



2vr 



+ OC 



h( T )e~ iuJT dr, 



(23) 
(24) 



the differential phase can be written as 
Acj)(t) = (1 + cos 9) cos 2a ■ 

/ + °° LO 
— ft,(w)sin(wZsin6»sin0)e iu;t dw. (25) 
-oo w 

The output of the detector, AcA(i), does not mirror 
the gravitational waveform, hit). Rather, the output is 
proportional to the convolution of the waveform with a 
square time-domain window whose amplitude and du- 
ration are determined by the orientation and baseline 
length of the interferometer. An alternative approach to 
the extraction of the gravitational wave signal from the 
interferometer's output is to differentiate the differential 
phase with respect to co-ordinate time. Adjusting the 
origin of the time co-ordinate in 122|) , and differentiating 
with respect to t, one finds 



dAcA 
dt 



(1 + cos 9) cos2a[h(t) — h(t — 2^sin#sin< 



Although this equation is in the more standard form for 
the description of the output of interferometric beam de- 
tectors !). equation (5)], this method of signal extraction 
makes noise estimation more complex and so this paper 
will confine itself to an examination of the phase differ- 
ence between the receiving stations alone. 



IV. THE RESPONSE TO SPECIFIC 
GRAVITATIONAL WAVEFORMS 

A. Monochromatic Gravitational Wave 

Here we shall consider the special case of a monochro- 
matic gravitational wave with angular frequency w gw , de- 
scribed in the frequency domain by 

h{w) = h- [S(ui — Wgw) + 5(uj + w gw )] . (26) 
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In the time domain this is equivalent to h(t) = h cos cu gw t. 
Substituting h(u) into l|25|l results in a time-averaged 
RMS differential phase of 



A(j) T 



(1 + cos 9) cos 2a sin(ui gw l sin 9 sin </>) | , 



V2 



(27a) 



where 



$ e ff = |(1 + cos 9) cos 2a sin(^ gw Z sin 9 sin <j>) \ . 

(. 27b ) 

<fr c g gives the differential phase per unit of gravitational 
wave strain and is a measure of the interferometer's sen- 
sitivity to gravitational waves. One can visualize $ g 
as follows. Consider replacing the interferometer with a 
simple distance measuring device consisting of two test 
masses and a (very precise) measuring stick marked off in 
radians of phase of the radio wavelength being observed 
by the interferometer. A passing gravitational wave will 
move the test masses back and forth, and using the mea- 
suring stick one can read off the change in the distance 
between them in radians. $ ff gives the number of ra- 
dians apart the test masses would need to be placed in 
order to observe the same change in their separation, in 
radians, as there is differential phase seen in the radio 
interferometer. $ e fj shall be referred to as the interfer- 
ometer's "effective phase length." The interferometer can 
be assigned a "natural" phase length of u> T 2l (the sepa- 
ration of the receiving stations in radians of the radio 
wavelength being observed), and then the ratio 



<l>r 



(28) 



which ranges from to e max ~ 0.65, can be thought of 
as the detector's efficiency. 

Since the efficiency is equivalent to the effective phase 
length normalized with respect to a combination of in- 
terferometer parameters, it is not useful for comparing 
the behaviour of different interferometers. Since the ef- 
ficiency is, in fact, only a function of the orientation of 
the gravitational wave with respect to the interferometer 
and the wavelength of the gravitational wave in units of 
the separation of the two receiving stations, it is a conve- 
nient way of discussing the dependence of any one inter- 
ferometer's sensitivity to the direction on the sky from 
which the gravitational wave is approaching, the orienta- 
tion of its polarization, and its wavelength. Plots of the 
efficiency showing its dependence on direction and wave- 
length are shown in Figure These show, specifically, 
the behaviour of the function 



f 



[ 1 + cos 9) cos 2a sin(w gw I sin 9 sin < 



both in the plane of the interferometer, = ±7r/2 (figures 
on the left), and over the full sky (figures on the right), 
for three choices of gravitational wave wavelength: from 
top to bottom co gw 2l = rad, 37r/2 rad, and Sir rad. In 
all plots, the polarization of the gravitational wave is cho- 
sen so as to maximize the efficiency (a = rad) . For the 
plots of the efficiency in the plane of the interferometer, 
the scale is such that the outer ring indicates an efficiency 
of I and rings are drawn at intervals of 0.1. The three 
filled circles indicate the orientation of the interferometer 
elements with the uppermost circle indicating the direc- 
tion toward the radio source (9 — rad). In the full-sky 
contour plots, the filled circle at the centre of each plot 
indicates the direction toward the radio source, while the 
outer ring corresponds to the direction away from the 
radio source and the two filled circles at ir/2 rad from 
the centre indicate the directions toward the receiving 
elements. Contours are drawn at —3 dB, —6 dB, and 
—9 dB relative to the detectors best possible efficiency of 
£max ~ 0.65. Several features are noteworthy: 

• The direction of greatest sensitivity to gravita- 
tional waves is not perpendicular to the interferom- 
eter's lines of sight to the radio source, where one 
might expect the wave's polarization to be most 
favourable, but rather lies in a direction displaced 
from the perpendicular toward the radio source. 

• The interferometer is almost completely insensitive 
to gravitational waves approaching from directions 
away from the radio source. 

• Compared to the response to longer wavelengths 
(Figure |2Ji), the interferometer responds poorly to 
gravitational waves with wavelengths shorter than 
the interferometer's baseline (Figure Et). 

To more clearly see how the frequency response of the 
interferometer depends on the baseline length, Figure 01 
shows the interferometer's effective length in the direc- 
tion of greatest sensitivity as a function of the gravita- 
tional wave's angular frequency for three choices of base- 
line. For the purpose of illustration, some fairly long 
baseline lengths were chosen: in terms of the light cross- 
ing time, they have lengths of 1 s, | s, and j s. The 
gravitational waves with these wavelengths have angular 
frequencies of 27rrad/s, 47r rad/s and 87rrad/s respec- 
tively. Recalling that the effective length is a measure 
of the sensitivity of the interferometer to gravitational 
waves, Figure 0] shows a flat response to gravitational 
waves with wavelengths longer than the interferometer's 
baseline. The poor short wavelength performance is at- 
tributable to the differential phase being the convolu- 
tion of the gravitational waveform with a square window 
whose duration is comparable to the light-crossing time 
of the interferometer's baseline: gravitational waves with 
periods less than this are smoothed out, diminishing their 
influence on the interferometer's output. Lest one think 
the solution to the poor high frequency performance is to 
use shorter baselines, we also see from Figure0]that going 
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(a) 



(b) 




FIG. 3: The efficiency of the interferometer's response to monochromatic gravitational waves at several wavelengths. From top: 
u)g V 2l = rad, 3tv/2 rad, and 3tv rad. The three filled circles in each plot indicate the orientation of the source and receiving 
elements. On the left are plots of the efficiency in the plane of the interferometer with the outer ring being an efficiency of 1 
and with rings drawn at intervals of 0.1. On the right are the same functions shown as contour plots over the full sky with 
contours drawn at —3 dB, —6 dB, and —9 dB relative to the detector's best possible efficiency of e max ~ 0.65 (visible as the 
peak efficiency in the top-leftmost plot). 
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FIG. 4: The interferometer's effective length (in the direction 
of greatest sensitivity) in units of light travel time as it de- 
pends on the angular frequency of the gravitational wave and 
the length of the interferometer's baseline. For the purpose 
of illustration, baseline lengths of 21 /c = 1 s (solid curve), ^ s 
(long dashes) and -r s (short dashes) were chosen. A larger 
effective length means the generation of a larger differential 
phase for the same gravitational wave amplitude and thus a 
more sensitive detector. 



to shorter baselines does not improve the high frequency 
response at all. It does flatten the response to higher 
frequencies but it accomplishes this by reducing the low 
frequency sensitivity, making it as poor as that at high 
frequencies. The conclusion is that longer baselines are 
better: they improve the low frequency response with- 
out, in fact, harming the high frequency response. The 
interferometer's flat response at low frequencies can be 
seen analytically by checking that 

lim -=^- = 0, 

w gw z^o dug W 

and the actual efficiency the interferometer asymptotes 
to at low frequencies is given by 

lim e = — |(1 + cos6>) sin sin ^ cos 2a | . 

LJg w / — >0 2 

It is interesting that although longer baselines do im- 
prove the overall sensitivity of the interferometer, con- 
trary to what one might expect upon consideration of 
Figure n it is not necessary to construct an interferome- 
ter with a long baseline in order to avoid being insensitive 
to long wavelength gravitational waves. One might have 
thought, for example, that a gravitational wave with a 
wavelength much longer than the separation of the re- 
ceiving elements would have the two lines of sight always 
sitting at approximately the same phase within the wave, 
thus making it difficult to detect its presence. However, 
the long wavelength also means the gravitational wave's 
effect on the radio signal is integrated over a longer seg- 
ment of the line of sight. The two effects balance each 
other giving the interferometer a fiat response at low fre- 
quencies. 



Although this section has looked specifically at an ex- 
actly monochromatic gravitational wave, the results can 
be reasonably applied to many anticipated astrophysical 
sources of gravitational waves. From 122|) . the relation- 
ship between the differential phase and the gravitational 
waveform only involves the integration of the waveform 
over an interval roughly equivalent to the light-crossing 
time of the interferometer's baseline. The results of this 
section continue to hold so long as the gravitational wave 
is approximately monochromatic over this interval. For 
radio receiving elements restricted to the surface of the 
Earth, this is an interval less than 40 ms. Even for space- 
based antennas separated by millions of kilometres the 
interval is on the order of 10 s. Over these time scales, 
many anticipated astrophysical sources of gravitational 
waves can be approximated as exactly periodic. 



B. Black Hole Ring-Down Gravitational Wave 

One example of an astrophysical system that can't rea- 
sonably be approximated as exactly periodic over tens of 
milliseconds is a perturbed black hole. When a spin- 
ning black hole is distorted, the irregularities of its event 
horizon are radiated away in the form of gravitational 
waves. The particular waveform generated by this pro- 
cess is called a "black hole ring-down" waveform and is 
of the form [HEl 



h(t) 



0, t < s, 

her u e t/(2Q) s i nL j s t, t>0s, 



(29) 



where 



uj s « (200 krad/ sec) 1 - 0.63(1 - a)™ 



M 



(30) 

is the centre frequency of the black hole's fundamental 
quadrupolar mode, 



Q « 2(1-0)" 



(31) 



is the "quality factor," M/Mq is the black hole's mass in 
solar mass units, and a is the black hole's dimensionlcss 
angular momentum (it's physical angular momentum is 
aMG/c 2 ). From 122(1 . the interferometer's differential 
phase, Acj)(t) , is related to the gravitational waveform by 
integration over an interval of time. Adjusting the origin 
of the time co-ordinate in (|22|l so that this interval just 
encounters the start of the ring-down waveform at t = 0, 
the differential phase produced by this type of wave is 



0J r 



A0(i) = hf(l 



for t > 0, where 

^ start 



cos 9) cos 2a 



e -cW(2Q) smu , T dr 



max(0, t — 21 sin 6 sin ( 
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FIG. 5: The ring-down waveform for a black hole of mass 
M = 50Mq and dimensionless angular momentum a = 0.3. 
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signal, as it is with monochromatic gravitational waves, 
so there is no clear way to define a quantity similar to 
the "effective phase length." Comparing to the analy- 
sis of monochromatic gravitational waves in the previous 
section, the RMS time average of the function plotted 
in Figure is algebraically equal to the efficiency times 
the baseline, e2l (hence the choice of milliseconds for the 
units of the vertical axis). The RMS average over infi- 
nite time, however, would give a result of 0, so for grav- 
itational waves localized in time this procedure is not a 
good way to assess the strength of the interferometer's 
output. 

The interferometer's response is plotted for two choices 
of the interferometer's baseline: 6000 km and 12000 km 
(light crossing times of 20 ms and 40 ms respectively). 
The shape of the interferometer's response as a function 
of time can be understood as follows. Prior to t = s, 
the gravitational wave has not yet encountered the in- 
terferometer and there is no differential phase seen. At 
t = s, the ring-down wave train wave crosses the line 
of sight of the first receiving clement, resulting in that 
station accumulating a phase change in its received ra- 
dio signal. There is then a constant differential phase 
between the stations until the wave train crosses the line 
of sight of the second receiving element (a little less than 
either 20 ms or 40 ms later) . At that time the second 
station accumulates the same change in the phase of its 
received radio signal and the differential phase between 
the stations returns to 0. 



FIG. 6: The interferometer's response to the ring-down wave- 
form for two choices of baseline: 6000 km (solid curve) and 
12000 km (dashes). These baselines have light-crossing times 
of 20 ms and 40 ms respectively. 



This evaluates to 



A<p(t) 



Q 



0Js 1 + 4Q : 



■(1 + cos 6*) cos 2a 



-w b t/(2Q) 



(2Q coslu s t + sin lv s t) 



tsta 



(32) 



A detailed exploration of the response of the radio in- 
terferometer to the full parameter space of these ring- 
down gravitational waveforms is outside the scope of this 
document. To illustrate the nature of the interferome- 
ter's response, we consider the single choice of param- 
eters consisting of M = 50Af Q , and a — 0.3. Such a 
black hole's fundamental quadrupolar mode has a centre 
frequency of lo s = 1736rad/s, and the ring-down wave- 
form has a quality factor of Q = 2.348. The ring-down 
waveform is shown in Figure [S] The interferometer's re- 
sponse to this type of gravitational wave is depicted in 
Figure which is a plot of 2A(f)(t) / (hui T ) as a function 
of time when the gravitational wave is approaching from 
the optimum direction with the optimum polarization. 
For these ring-down waveforms, the output of the inter- 
ferometer is not proportional to the gravitational wave 



V. NOISE IN THE INTERFEROMETER 

To assess the viability of the type of gravitational wave 
detector described in this paper, it is necessary to investi- 
gate the noise that must be anticipated in the detector's 
output. A full analysis of the noise level requires a com- 
plete specification not only of the detection apparatus 
but also of the observation strategies and data analysis 
techniques to be used in the extraction of a gravitational 
wave signal from the detector's raw output. 

The technique of gravitational wave detection de- 
scribed in this paper can, in principle, be implemented 
with both optical and radio technologies. Although it 
has been implied throughout the document that it is ra- 
dio systems that are being considered, the analysis up to 
this point is applicable to both. Optical systems benefit 
from operation at higher angular frequencies but radio 
systems, at least terrestrially, benefit from much longer 
baselines and the ability to record the electromagnetic 
signals for detailed off-line analysis. Another advantage 
of the implementation of the detector using radio tech- 
nology is the possibility of simultaneous multi-use of the 
technology for other scientific research. We now specif- 
ically assume a radio-based system. There will be little 
discussion of observation strategies or data analysis tech- 
niques as these are beyond the scope of this paper. 

Gravitational wave detection proves to be one of 



11 



Im(signal) 
signal. 




noise 

signal + noise 



Re (signal) 



FIG. 7: Complex signal and noise phasors for the radio signal 
being received by the interferometer. 



the most challenging measurement problems in modern 
physics. All detection techniques currently being used or 
considered require extraordinary technological achieve- 
ments. The technique described here does not provide 
a "magic bullet" for circumventing any of the difficul- 
ties in gravitational wave detection so one must expect 
that the practical implementation of this technique will 
also require an apparatus with extraordinary operating 
parameters. 



A. Thermal Noise in the Electronics 

The antenna structures, radio receivers and data ac- 
quisition systems used in astronomical radio intcrfcrom- 
etry applications are not noise-free, and contribute to 
the uncertainty in the measurement of the interferometer 
phase. Figure shows an example complex phasor dia- 
gram for the radio signal being received and the thermal 
noise added to it by the interferometer. If the signal-to- 
noise ratio is greater than 1, then it can be seen that the 
phase error introduced by the noise is 



noise 



| signal 



= SNR7 1 



where the subscript r on the right-hand side identifies 
this as the radio interferometer's signal-to-noise ratio. If 
we take our gravitational wave signal to be the time- 
averaged RMS differential phase in (|27a|l and compare 
this to the phase noise from above, we get 



SNR„ 



A0 r 



A</> rros SNR r 



as the signal-to-noise ratio on a single determination of 
the differential phase due to the presence of a gravi- 
tational wave. Since Acf> rms oc hrjf 1 - is typically much 
smaller than 1, we see immediately that we are going to 
need a very high signal-to-noise ratio in the interferome- 
ter in order to see gravitational waves above the thermal 
noise. 



The signal-to-noise ratio of a two-element radio inter- 
ferometer is given by 0, page 158] 



SNR r 



/2tt 



' TaiTa 2 

2sy sl T s y S2 



where Aw r / (2tt) is the bandwidth of the receiving system 
in Hertz, At is the coherent integration time in seconds, 
TXi is the antenna (signal) temperature in Kelvin at the 
i th antenna, and T sySi is the system (noise) temperature 
in Kelvin at the i th antenna. The antenna temperature 
is given by [l3| 



T A = 



10- 26 ^A cff 
2fc R 



where S v is the spectral flux density of the observed ra- 
dio source in Janskys (1 Jy = 10~ 26 Wiir 2 Hz _1 ), A c s 
is the antenna's effective aperture area in square metres, 
and fee = 1-4 x 10~ 23 JK _1 is the Boltzmann constant. 
The factor of 2 is due to the radio receiver being able 
to observe only one electromagnetic polarization. For a 
"blind search" observation of gravitational waves with an 
angular frequency of w gw , our interferometer's coherent 
integration time cannot exceed the inverse Nyquist fre- 
quency of At = 7r/ti> gw . If we, further, assume the two 
receiving stations are identical (same effective apertures, 
same system temperatures), then altogether the phase 
uncertainty due to thermal noise in the radio interferom- 
eter's electronics is 



4x 10 3 



T 



Aw r S u A cS 



(33) 



where 5' v is in Jy, T S y S is in K, and ^4eff is in m 2 



B. Other Noise Sources 

The thermal noise in the interferometer's electronics 
limits one's ability to measure with precision the dif- 
ference in phase between the signals being received at 
each station. Even given an exact measurement of the 
interferometer's differential phase, it is still difficult to 
attribute it exclusively to a passing gravitational wave as 
there are other processes that contribute to phase differ- 
ences between the receiving stations. Examples of such 
systematic noise sources are: 

• Time-dependant deformations of the interferometer 
due to, for example, seismic noise and tidal distor- 
tion of the Earth. 

• Variations in the refractive index of the medium 
through which the radio waves travel on their way 
to the receiving stations. These include variations 
in the refractive index of the Earth's atmosphere as 
well as possibly interplanetary or interstellar gases 
(depending on the distance to the radio source). 
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• Thermal distortions of the antennas and the conse- 
quent displacement of the antenna phase centres. 



• Variations in the instrumental signal delays 
through cables and electronics at each antenna site. 



These noise sources, as well as others not mentioned 
here, affect present day radio astronomical and geode- 
tic radio interferometric observations. In general it will 
be necessary to either remove the effects of systematic 
noise sources by clever observing strategies and/or so- 
phisticated signal processing techniques, or to provide 
calibration systems that are capable of accounting for 
their presence. In some cases, present-day radio interfer- 
ometric techniques may be capable of providing the nec- 
essary tools. For example, dispersive effects of the iono- 
sphere and the interplanetary/interstellar medium can 
be removed by multi-frequency observations and radio 
interferometry phase calibration (peal) systems fl2l page 
290] can remove the effects of variations in instrumental 
phase delay between stations at the picosecond level ^ij . 

Differential observing strategies may possibly be devel- 
oped that are effective in removing systematic effects in 
the output data. Figure |3| shows that, as a gravitational 
wave detector, a radio interferometer is reasonably direc- 
tional, and the use of antenna arrays and beam-forming 
techniques might be employed to implement a gravita- 
tional wave detector with high directivity. Such a sys- 
tem combined with multi-beam radio frequency anten- 
nas could, in principle, carry out a differential search for 
gravitational waves incident from several directions on 
the sky simultaneously. Differential observing strategies 
such as this are capable of cancelling systematic noise 
sources to a high precision. 

The sources of noise mentioned above produce differ- 
ential phase in the interferometer through means other 
than the curvature of space. Even when all such noise 
can be removed from the interferometer's output, leav- 
ing only curvature-induced differential phase, there are 
sources of curvature other than gravitational waves. Phe- 
nomena such as seismic and atmospheric activity can pro- 
duce changes in local gravity comparable to the expected 
amplitude of astrophysical gravitational waves @. Unlike 
ground-based laser interferometers, the radio signals used 
by the type of detector investigated here do not travel 
along geodesies that are entirely local to the Earth. For 
this reason, it is not yet clear the extent to which such 
sources of noise contribute to the detector's output but 
it is possible this type of detector is less sensitive to them 
than are entirely ground-based systems. 

A full analysis of the noise from these and other sources 
and the extent to which these can be subtracted from the 
detector's output will be the subject of future work. 



VI. PRACTICAL GRAVITATIONAL WAVE 
DETECTION 

If we assume for the moment that all systematic 
sources of differential phase can be exactly accounted 
for and subtracted from the detector's output, then that 
leaves the interferometer's thermal noise as the only im- 
pediment to the operation of this type of detector. While 
this is an unrealistic situation, it does represent a sensi- 
tivity limit which cannot be improved upon through any 
amount of creative observation strategies or signal ex- 
traction techniques (with the exception of extending the 
interferometer's integration time, and this will be dis- 
cussed below). 

For good gravitational wave signal to noise, it has been 
established that long baselines, high frequency radio ob- 
servations, low receiver noise temperatures and bright 
radio sources are required. In addition, one can con- 
sider performing multiple simultaneous observations of 
a gravitational wave: by constructing an interferome- 
ter with more than two elements, the differential phase 
across multiple baselines can be measured. If these mea- 
surements are time-shifted appropriately and summed to- 
gether, the thermal noise in the individual measurements 
will sum incoherently while any signal that is present 
will be summed coherently. For N baselines, the relative 
phase noise could be reduced by as much as a factor of 

s/n. 

Setting A0 rms in l|27a|l equal to 8cj) x in (|33[1 and solv- 
ing for h gives us the gravitational wave shear needed to 
achieve a signal-to- noise ratio of at least 1. This is 

h _fiv 111 3 / ^gw Tsys 

Vin "OX 1U W — — — — • 

y U^AWr b„A e g 

1(1 + cos 9) cos 2a sin(w gw i sin 9 sin <p)\~~ . (34) 

Optimistic yet achievable numbers for some of the pa- 
rameters in this expression are: a radio observation cen- 
tre frequency of u> T = 10 11 rad/s, a bandwidth of Aw r = 
10 9 rad/s, system temperatures of T sys = 20 K, effective 
apertures of A e g = 10 3 m 2 , and an interferometer con- 
sisting of 32 stations for a total of 496 baselines of about 
21 = 2000 km (1 — 3 ms) each. Given these parameters, 
Figure |H1 shows plots of h nlin for radio source spectral flux 
densities of 10 2 Jy and 10 4 Jy (at 16 GHz). Flux densi- 
ties of this order are available naturally from galactic wa- 
ter vapour maser sources at 22 GHz or could be provided 
artificially, for example by means of a spacecraft-born ra- 
dio transmitter. The Canadian Large Adaptive Reflector 
(CLAR), embodies a design for a low-cost, large aperture 
(300 m diameter or larger), fully-steerable, multi-beam 
radio astronomy antenna fl7l fl8| . One technical im- 
plementation for the Square Kilometre Array (SKA) en- 
visages the construction of an array of CLAR antennas. 
Such an interferometer array might, potentially, have use 
in the search for gravitational radiation although its an- 
ticipated operating parameters differ somewhat from the 
numbers quoted above. 
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FIG. 8: The minimum shear required in order to achieve a 
signal to noise ratio of 1 in the observation of gravitational 
waves using a radio interferometer array with the parameters 
listed in the text. Two radio source spectral flux densities 
are considered: 10 2 Jy (solid line) and 10 4 Jy (dashed line). 
Compare to Figure 4 of [T^ . 

To compare the sensitivity of this type of detector to 
those of other types, as well as to the strengths of antic- 
ipated sources of g ravitational waves, see the recent re- 
views found in |l5l | and The noise curves in Figure 
are well above any anticipated gravitational wave signals. 
It must be recalled, however, that these noise curves are 
for the raw detection of a gravitational wave signal at or 
above the noise level: a pen on a chart recorder marking 
out the gravitational wave shear as a function of time. 
This picture is not what is meant when it is said that 
the current generation of gravitational wave observato- 
ries will "observe" gravitational waves. Rather than ob- 
taining a direct read-out of the gravitational wave shear 
as a function of time, one convolves the output of the 
detector against a template waveform and compares the 
correlation coefficient to a baseline coefficient obtained 
by correlating the same template against a model of the 
detector's noise. By this method, one extends the coher- 



ent integration time of the interferometer and reduces 
the output to the probability that the given waveform is 
present in the data. In the case of many astrophysical 
sources, such as binary neutron star systems, integration 
times of months or even years are anticipated |19| . This is 
to be contrasted to the integration times of < 10 3 s that 
were assumed for the computation of Figure [SJ Since the 
noise curves in Figure [S] are inversely proportional to the 
square root of the integration time, then if such long in- 
tegration times could be used in the radio interferometer 
one could potentially decrease /i m in by two or three orders 
of magnitude. In the frequency band shown in Figure |S1 
such a reduction in the noise levels would make this type 
of detector competitive with the LISA detector's noise 
level as shown in Figure 4 of [l5| . 

It must be stressed, however, that systematic noise has 
been ignored. An analysis of the contribution from sys- 
tematic noise is complicated by the availability of sophis- 
ticated calibration techniques that can subtract a great 
deal of it. A complete analysis of the influence of sys- 
tematic noise will involve not only an estimate of the raw 
contribution from these sources but also the development 
and subsequent analysis of the effectiveness of observa- 
tion strategies and data processing techniques. This will 
be the subject of future work. 



VII. CONCLUSIONS 

The technique, introduced here, of using the time- 
dependent phase difference between two radio receiving 
stations to watch for the passage of a gravitational wave 
would seem to be a promising method of gravitational 
wave observation worthy of further investigation. In the 
future, a more detailed analysis of realistic interferome- 
ter configurations and data analysis techniques will shed 
more light on true sensitivity of this type of detector. 
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